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ABSTRACT 



The low mass X-ray binary (LMXB) associated with the 
M31 globular cluster Bo 158 is known to exhibit inten- 
sity dips on a —2.78 hr period. This is due to obscura- 
tion of the X-ray source on the orbital period by material 
on the outer edge of the accretion disc. However, the 
depth of dipping varied from <10% to —83% in three 
archival XMM-Newton observations of Bo 158. Previ- 
ous work suggested that the dip depth was antic orrelated 
with the X-ray luminosity. However, we present results 
from three new XMM-Newton observations that suggest 
that the evolution of dipping is instead due to precession 
of the accretion disc. Such precession is expected in neu- 
tron star LMXBs with mass ratios <0.3 (i.e. with or- 
bital periods <4 hr), such as the Galactic dipping LMXB 
4U 1916— 053. We simulated the accretion disc of Bo 
158 using cutting-edge 3D smoothed particle hydrody- 
namics (SPH), and using the observed parameters. Our 
results show disc variability on two time-scales. The disc 
precesses in a prograde direction on a period of 8 1 ±3 hr. 
Also, a radiatively-driven disc warp is present in the inner 
disc, which undergoes retrograde precesson on a —3 1 hr 
period. From the system geometry, we conclude that the 
dipping evolution is driven by the disc precession. Hence 
we predict that the dipping behaviour repeats on a —81 hr 
cycle. 

Key words: X-rays: general; X-rays: binaries; Galaxies: 
individual: M31, Accretion: accretion discs. 



iTrudolvubov et al]j2002l) report —83% modulation in the 
0.3-10 keV flux of Bo 158 on a 2.78 hour period dur- 
ing the —60 ks 2002 January XMM-Newton observation. 
The modulation resembles the intensity dips seen in high 
inclination LMXBs due to photo-electric absorption of 
X-rays by mat erial that is raised ab ove the body of the 
accretion disc (White & Swankl H982l) . They also report 
—30% dips in the 2000, June XMM-Newton lightcurve 
and -50% dips in the 0.2-2.0 keV lightcurve of the 1991, 
June 26 ROSAT/PSPC observation. However, no sig- 
nificant dips were found in the 0.3-10 keV lightcurve 
of the 2001 June XMM-Newton observation; the au- 
thors placed a 2a upper l imit of 10% on the modulation. 
Trudolvub ov et al.l 112002) concluded that the depth of the 
intensity modulation was anti-correlated with the source 
luminosity. 

We present further XMM-Newton observations and mod- 
elling results which suggest that the variation in dipping 
behaviour may instead be due to precession of the accre- 
tion disc. Such behaviour is associated with the "super- 
hump" phenomenon that is observed in interacting bina- 
ries where the mass r atio of the secondary to the primary 
is smaller than —0.3 Jwh itehurst & Kingtri991l) . Super- 
humps are briefly reviewed in Sect. [2] followed by details 
of the observations and data analysis in Sect. [3] and our 
results in Sect. [4] Numerical modelling of the system 
is discussed in Sect. 5; the system was simulated by a 
3D smoothed particle hydrodynamics (SPH) code. We 
present our discussion in Sect. [6j and finally our conclu- 
sion in Sect. 7. 



1. INTRODUCTION 



2. SUPERHUMPS 



Bo 158 is source number 158 in the cata logue of globu- 
lar clu sters that were identified in M31 by Battistini et al. 
( 1987). Its X-ray counterp art is located at a = 
QQ h 43 m 14.2 s , S = 41°07'26.3" tDi Stefano etaill2002l) . 
Trudolvu bov et alJ ( 120021) identified the X-ray source as 
a likely low mass X-ray binary (LMXB) with a neutron 
star primary; following their work, we will use the desig- 
nation "Bo 158" to describe the X-ray source here. 



Superhumps were first identified in the superoutbursts of 
the SU UMa sub-class of cataclysmic variables. They are 
manifested as a periodic increase in the optical brightness 
on a period that is a few pe rcent longer than the orbital pe- 
rio d JVogtlll974tlWarnerl[T975T) . In the model proposed 
by Osaki ( 1989), superhumps occur when the outer disc 
reaches a 3:1 resonance with the secondary. The addi- 
tional tidal forces exerted on the disc by the secondary 
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Table 1. Journal of XMM-Newton observations of the 
M31 core. A1-A3 are available in the public archive, 
while P4-P6 are proprietary observations. Two observa- 
tions were made during 2004 July 19; the first observa- 
tion (a) started at 01:42:12, and the second (b) started at 
13:11:22. 



Observation 


Date 


Exp 


Filter 


Al 


2000 Jun 25 


34 ks 


Medium 


A2 


2001 Jun 29 


56 ks 


Medium 


A3 


2002 Jan 6 


61 ks 


Thin 


P4 


2004 Jul 17 


18 ks 


Medium 


P5 


2004 Jul 19a 


22 ks 


Medium 


P6 


2004 Jul 19b 


27 ks 


Medium 



at this stage cause the disc to elongate and precess, and 
also greatly enhance the loss of angular momentum, in- 
creasing the the mass-transfer rate. The disc precession is 
prograde in the rest frame, and the secondary repeats its 
motion with respect to the disc on the beat period between 
the orbital and precession periods, slightly longer than the 
orbital period. The secondary modulates the disc's vis- 
cous dissipation on this period, giving rise to maxima in 
the optical lightcurve, known as superhumps. 

The requirement for the 3:1 resonance to fall within the 
disc's tidal radius is that the mass r atio of the secondary 
to the primary be less than ^0.33 (Whitehurst & King, 
II991T) . If we assume that the secondary is a main se- 
quence star that fills its Roche lobe, then the relation mi 
~ 0. 1 1 Phr holds, where to 2 is the mass of the secondary 
in solar units and Ph r is the orbital period in hours (e.g. 
iFrank. King. & Raind 120021) . Hence any accreting bi- 
nary that has a short enough orbital period may exhibit 
superhumps. 



3. OBSERVATIONS AND DATA ANALYSIS 



In addition to the three XMM-Newton observations anal- 
ysed by iTrudolvubov et alJ (12002), we conducted a pro- 
gramme of four ^20 ks observations over 2004, July 16- 
19. We present results from our analysis of the archival 
data, along with three of the four 2004 observations; 
the other observation suffered from flaring in the parti- 
cle background over 90% of the observation and is not 
considered further here. A journal of observations is pre- 
sented in Tabled 

We analysed data from the pn, MOS 1 and MOS2 instru- 
ments, which share the same 30" x 30" field of view. For 
each observation, we selected a circular extraction region 
with a 40" radius, centred on Bo 158, and an equivalent 
source-free region for the background. The background 
region was on the same chip as the source, and at a sim- 
ilar angular offset from the optical axis. We extracted 
lightcurves from the source and background regions in 
the 0.3-10 keV energy band with 2.6 s binning, and also 
obtained pn spectra of the source and background regions 



along with the associated response files. The spectra were 
then grouped for a minimum of 50 counts per bin. 



4. RESULTS 



The 0.3-10 keV EPIC (MOS1 + MOS2 + pn) lightcurves 
of Observations A1-A3 are presented in Fig.^ with x- 
and y- axes set to the same scale, and with 200 s binning. 
Most striking is Observation A3, wit h six dipping inter- 
vals o n a 10017±50 second period (ITrudolvubov et all 
120021) : the structure of the dipping is seen to vary sub- 
stantially even over one observation o f 60 ks. We e mpha- 
sise th at the 30% dipping reported by Trudolvub ov et alJ 
(2002) for observation Al is an overall average; a deep 
dip is seen at ^28 ks into the observation, but very lit- 
tle evidence of dipping is observed in the intervals of 
expected dipping at ^8 ks and ^18 ks into the observa- 
tion. In this regard, Observation Al resembles the 1985, 
October lightcurve of the Galactic superhumping LMXB 
4U 1916-053, observed by EXOSAT. In that observa- 
tion, de ep dips are observe d only after the first four orbital 
cycles (S male et all 11988b . Little evidence of variability 
is seen in the lightcurve of Observation A2. 

In Fig. [21 we present the 0.3-10 keV EPIC lightcurves of 
Observations P4-P6; the x-axis is scaled to P6, and the y- 
axis matches that of Fig.^ The most prominent feature 
is the dip in P5; it has a depth of ^100% and a total dura- 
tion of ^2500 s. Using the period of Trudolvub ov et alJ 
(2002), we identified the expected times of dipping, la- 
belled 'D', in P4 and P6, using the deepest part of the dip 
in P5 as phase zero. 

In the P4 lightcurve, there is no evidence for dipping 
during the first expected dip interval, but some evidence 
of dipping ~4000 s after the second interval, 20 cycles 
away from the dip in P5. Hence, this possible dip would 
require a period that is either ~200 s shorter or ~320 
s long er than the 10017 s given by ITrudolvubov et alJ 
(2002). However, there is no other evidence for these 
other periods in the lightcurves of P4, P5 or P6; hence it 
is clear that the dipping behaviour of Bo 158 evolves on 
a time scale of a few days. 

Several spectral models were fitted to the 0.3-10 keV pn 
spectrum of P4, each suffering absorption by material in 
the line of sight. P4 was chosen because it had the longest 
interval of persistent emission that was not contaminated 
by background flares; the resulting source spectrum con- 
tained ~2200 cou n ts. W e applied the two models that 
ITrudolvubov et all ( 120021) used to model the spectra of 
A1-A3, namely a power law model and a Comptoni- 
sation model (COMPTT in XSPEC). We also applied a 
two component model. The emission of many Galac- 
tic LMXBs has been succesfully described by a model 
consisting of a blackbody and a cut-off power law (e.g . 
Church & Balucihska-Church, 1995: Church et all ! 1998: 
Barna rdet all |2003); we approximate this model to a 
blackbody + power law model, because of the narrow 
pass-band. 
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Figure 1. Combined EPIC 0.3-10 keV lightcurves (in 
count s" 1 ) of Bo 158 from the archival XMM-Newton 
observations, A1-A3; the dates of each observation are 
shown. 



We find that the best fit parameters for the power law and 
COMPTT models agree w ell with the values presented by 
iTrudolvubov et alJ d2002h : however, the two-component 
model provided the best fit. We find a 0.3-10 keV flux 
of ~ 2xl0 -12 erg cm~ 2 s _1 for all the fits to the P4 
data; this gives a 0.3-10 keV luminosit y of ~1.4xl0 38 
erg s~ 1 , assuming a distance of 760 kpc llvan den Bergh. 
2000). 



The depth of dipping in Al varies from ~0 to ^70% with 
no significant change in the mean intensity, suggesting 
that the amplitude of dipping is not simply anticorrelated 
with the source luminosity. Instead, the variation in dip- 
ping behaviour may be caused by disc precession. This 
hypothesis motivated our simulation of the accretion disc 
in Bo 158, using three dimensional smoothed partical hy- 
drodynamics, discussed in Sect. [5] 



The lightcurves of XMM-Newton observations of Bo 158 
show no true eclipses; hence, we know that we are not 
viewing the system edge on. If the disc were tilted with 
respect to the binary plane, and precessing, then one 
might expect to observe dips in some part of the disc 
precession cycle, but not in others. Dipping is observed 
throughout observation A3; this suggests that the dipping 
phase in the disc precession cycle lasts > ~60 ks. The Al 
lightcurve covered three intervals of expected dipping, 
yet only one dip is seen, toward the end of the obser- 
vation; we suggest that this dip signals the onset of the 
dipping phase. Contrariwise, P5 and P6 appear to sample 
the end of the dipping phase, as a dip is observed in P5, 
yet no dips are seen in P6. 



.1111111 

, , , , i , , 


■ * i * 

. . i . 




2004 Jul 

i i i i i i i 


1 l 1 ■ 
17 : 

ill' 


. . . . . i 

t * 




. . i . . . . 


2004 Jul 


19aj 




1 


4 




2004 Jul 

M 





5000 10 4 1.5xl0 4 2xl0 4 2.5xl0 4 

Time (s) 



Figure 2. Combined EPIC 0.3-10 keV lightcurves (in 
count s~ l ) of Bo 158 from the 2004 XMM-Newton obser- 
vations, P4-P6; the dates of each observation are shown. 
The horizontal line in the P5 lightcurve indicated the in- 
terval used for spectral fitting. The shaded area of the 
P6 lightcurve indicates a period of background flaring. 
Times when dipping is expected are labelled "D", using 
the deep est part of the dip in P5 as time zero, and the 
period o nTrudohubov et al\ §2002). 

5. SPH SIMULATION OF THE DISC 



5.1. Binary Parameters 

The accretion disc was modelled using a three- 
dimensional Smoothed Particle Hydrodynamics (SPH) 

co mputer c ode tha t has been describ e d in detail 

in iMurravl dl996l 119981) iTruss et alJ d2000l) . and 
Foulkes, Haswell, & Murray ( 200gh We assumed 
the orbital period to b e 10017 s, as obtained by 
ITrudolvubov et alJ J2002I) . 

The dipping source Bo 158 is a bright globular cluster X- 
ray source, with a 2.78 hr binary period. Thirteen Galac- 
tic globular clusters contain bright X-ray sources; twelve 
of these are neutron star LMXBs, whi le the primary o f 
the other one is unknown (see e.g. lin't Zand et all 120041) . 
Hence Bo 158 is a likely neutron star LMXB, and we as- 
sume the primary mass to be 1.4 M Q . 

For the secondary, we considered a main sequence star 
and a white dwarf, sin ce 4U 1916— 053 has a likely white 
dwarf secondary (e.g. IChou et all 1200 ll) . For a Roche 
lobe-filling main sequence star the app roximate relation 
m ~ 0.11 Phr holds (e.g. iFrank etail 120021) . giving a 
mass of ^0.30 M Q . If instead th e star is a white dwarf, 
using the mass radius relation of Nauenberg ( 1972) and 
Roche geometry gives an implausibly small mass of 
0.005 Mq. Assuming a main sequence secondary, we 
found the mass ratio to be 0.2, indicating that superhumps 
and disc precession were likely. 

Finally, the luminosity of the system was taken to be 
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1.4xl0 38 erg s" 1 , the 0.3-10 keV luminosity of Bo 158 
in Observation P4. Such a high luminosity may be ex- 
pected to cause warping of the accretion dis c, even for a 
previously flat disc (see e.g. Pringle, 1996); warping is 
discussed in Sect. 5.3. 

The accretion disc had an open inner boundary condition 
in the form of a hole of radius r\ — 0.025a, where a is 
the binary separation, centred on the position of the pri- 
mary object. Particles entering the hole were removed 
from the simulation. Particles that re-entered the sec- 
ondary Roche lobe were also removed from the simula- 
tion as were particles that were ejected from the disc at a 
distance > 0.9a from the centre of mass. 

We assumed an isothermal equation of state and that 
the dissipated energy was radiated from the point at 
which it was generated, as e lectromagnetic radiation. The 
IShakura & Sunvaevl dl973l) viscosity parameters were set 
to ai ow = 0.1 and ahigh = 1.0, an d the viscos it y stat e 
changed smoothly as described in iTruss et all J2000I) . 
The SPH smoothing length, h, was allowed to vary in 
both space and time and had a maximum value of 0.01a. 



5.1.1. The gas stream 

We simulated the mass loss from the secondary by intro- 
ducing particles at the inner Lagrangian point (Xi). The 
mass transfer rate and the particle transfer rate were pro- 
vided as input parameters, and the mass of each particle 
was derived from these parameters. A particle was in- 
serted with an initial velocity (in the orbital plane) equal 
to the local sound speed of the donor, cd, in a direction 
prograde of the binary axis. The z velocity of the inserted 
particle was chosen from a Gaussian distribution, with a 
zero mean and a variance of O.Icd- However, the infla- 
tion of the site of collision between the gas stream and 
outer disc was not modelled. 



5. 1 .2. The initial non-warped accretion disc 

The simulation was started with zero mass in the accre- 
tion disc and with the central radiation source switched 
off. A single particle was injected into the simulation 
every O.OIQ^A at the L\ point as described above until 
a quasi-steady mass equilibrium was reached within the 
disc. This was taken to be when the number of particles 
inserted at the L\ point, the mass transfer rate, was ap- 
proximately equal to number of particles leaving the sim- 
ulation at the accretor, the accretion rate. The simulations 
were continued for another 3 orbital periods to ensure 
mass equilibrium. The number of particles in the sim- 
ulated accretion disc was approximately 40,000 giving a 
good spatial resolution; the average number of 'neigh- 
bours', i.e. the average number of particles used in the 
SPH update equations, was 8.2 particles. The simulated 
disc encountered the Lindblad 3:1 resonance and became 



eccentric. The disc precessed in a prograde direction giv- 
ing rise to superhum ps in t he si mulated dissipation light- 
curves (c.f. iFoulkes et all 120041) . The radiation source 
was then turned on which gave rise to a very small num- 
ber of particles being ejected from the accretion disc. 



5.2. Surface finding algorithm & self-shadowing 

Accretion-powered radiation from the inner regions of 
the disc and the accreting object itse lf exert a force on 
the irradiated disc surface. Following Pringle ( 1996) the 
radiation source is modelled as a point source at the cen- 
tre of mass of the accretor. To apply this force, particles 
on the surface of the accretion disc had to be identified. 
We used a convex hull algorith m to fi nd th e surface par- 
ticles as described in iMurravl d!998ft and IFoulkes et alJ 
(2006). A ray-tracing algorithm was used to determine 
regions of self-shadow. For each particle found on the 
disc surface a light-ray was projected from the particle to 
the position of the radiation source at the centre of the 
disc. The particle was deemed to be illuminated by the 
radiation source if this light-ray did not intersect any disc 
material between the particle surface position and the ra- 
diation source (i.e. the particle could see the central ra- 
diation source). The radiation force was only applied to 
particles that were considered to form part of the disc sur- 
face and were illuminated by the central radiation source. 



5.3. Disc warping and precession measure 

For an optically thick disc, a warp can develop as 
a result of the radiation force (e.g. iPringld [1 996; 
iQgilvie & Dubusll2001tlFoulkes et alll2006l) . This is due 
to the fact that any radiation absorbed by a specific region 
on the disc surface will be later re-radiated from the same 
spot, normal to the disc surface. Hence any anisotropy 
in the disc structure will cause an uneven distribution of 
back-reaction forces on the disc surface, further perturb- 
ing the disc. A sufficiently high luminosity can induce 
and su stain a warp even in an originally flat disc ( Pringle, 
I1996L and references therein). 



5.4. Numerical modelling results 

As a result of the disc precession, viscous stresses in the 
disc vary significantly with time. The resultant energy 
dissipation in different regions of the disc, and hence the 
disc luminosity, varies with time. The disc luminosity 
was not modelled in detail. We assumed that the lumi- 
nosity was directly related to the disc regions with sig- 
nificant energy release through viscous dissipation. The 
viscous dissipation heats the gas in the accretion disc and 
it is assumed that the heat is radiated away from the point 
at which it was generated. Superposed on a steady signal 
there is a repeating series of "humps." The spacing of the 
humps corresponds to the superhump period, but are not 
representative of true optical lightcurves. 
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Figure 3. Particle projection plots for the SPH model. The position of each particle is indicated by a small black dot. 
The left panel is a plan view of the accretion disc as seen from above the disc. The cross at the centre of the plot shows 
the position of the primary object. The solid dark line is the Roche lobe of the primary and the L \ point is to the right and 
middle of the plot. The right panels are are particle projection plots on a plane perpendicular to the orbital plane and 
through the system axis in the xz and yz directions. 



In order to determine the superhump period, P s h, we ob- 
tained a power density spectrum from ~30 superhump 
cycles of the simulated light curve. We estimated the 
superhump period to be (1.035 ± 0.005)P or b- This im- 
plies the precession period of the outer regions, P prcc = 
(29±l)P orb , or81±3hr. 

Fig. 13 shows a snapshot from the simulation, from sev- 
eral different angles. The left panel shows an x-y pro- 
jection of the disc. The disc is clearly asymmetric and 
elongated, and a spiral density wave is clearly seen. This 
wave is so intense that it is removing material from the 
accretion disc and return ing it back to the Roche lobe of 
the secondary; see Foulke s et alJ (120041) for a full detailed 
description of a similar system with a mass ratio of 0.1. 

The two upper right-hand plots of Fig. |5J labelled yz- 
view and xz-view, are side views of the disc in the y-z 
and x-z directions respectively. The disc warp is clearly 
apparent in these two plots. The warp is odd symmetrical 
about the centre of the disc. 

The warp amplitude and size precessed as a solid body 
in a retrograde direction relative to the inertial frame. We 
found that P warp ~ 1 1 P or b- Fig.|4]shows the radial profile 
of the warp for five consecutive orbital cycles; the maxi- 
mum extent of the warp is 6-11° above the plane of the 
disc. 



6. DISCUSSION 



Disc precession is inferred from the 0.3-10 keV 
lightcurves of Bo 158, as is expected given its extreme 
mass ratio (short orbital period). As such, it resembles 
the Galactic superhumping LMXB 4U 1916-053. Since 



the LMXB Bo 158 is in a globular cluster near the centre 
of M31, it is unlikely that the optical period will ever be 
known. However, our Fourier analysis of the simulated 
dissipation lightcurves indicates a superhump period that 
is 3. 5 ±0.5% longer than the orbital period. Given the 
association between the d ips and superhump period re- 
ported bv lRetter et alJ l l2002t) . the 10017 s period may be 
the superhump period, in which case, the orbital period 
would be ^4% shorter. Such shortening of the period 
would not dramatically affect the outcome of our SPH 
modelling. 



Our simulations of the disc show two distinct types of 
variability in the disc structure. First is the elongation 
and prograde precession of the disc due to tidal interac- 
tions with the secondary at the 3:1 resonance; the disc 
precesses on period of 81 ±3 hr. We also see warping of 
the accretion disc, driven by irradiation of the disc surface 
by the central X-ray source; the warp is stable exhibits 
retrograde precession on a ^3 1-hr period. 



It is therefore important to establish which region is re- 
sponsible for the observed variation in dip morphology. 
The lightcurves of observations A1-P6 show no eclipses. 
From Kepler's law and t he ratio of the secondary ratio 
to the binary separation (Easleton, 1983]), the secondary 
has an angular radius of ~15°; hence the inclination > 
^75°. We see from Fig. 6 that the disc warp does not de- 
viate from the plane of the disc by more than 1 1° in our 
simulations, suggesting that the observed dips are likely 
to evolve on the disc precession period. 
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Warp measure for five orbital periods 
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Figure 4. Radial warp profiles. The vertical axis is the 
warp amplitude. The horizontal axis is distance from the 
primary object normalized such that the binary separa- 
tion is 1. Plots (a), (b), (c), (d) and (e) are for five con- 
secutive orbital cycles. 

7. CONCLUSIONS 



We have analysed three new XMM-Newton observa- 
tions of the M31 dipping LMXB Bo 158, in addi- 
tio n to re-analysing the th ree observations discussed 
in iTru dolvubov et all d2002l) . The newer observations 
spanned ~3 days in 2004, July. We find that that the rela- 
tionship between source intens ity and depth of dipping i s 
not so simple as described by Trudolvu bov et alf (2002). 
Instead, we believe that the observed variation in dipping 
behaviour is caused by precession in the accretion disc; 
dipping would be confined to a limited phase range in the 
disc precession cycle. 

We modelled the accretion disc with 3D SPH, and found 
prograde disc precession on a 81 ±3 hr period, as well 
as radiatively driven disc warp that precessed on a 31 
hr period in a retrograde fashion. We find that the disc 
precession is most likely to affect the observed dipping 
behaviour. Hence, we predict that the dipping behaviour 
of Bo 158 experiences a 81 ±3 hour cycle; this period is 
consistent with the observed variation of the dips. 
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